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Abstract 

We analyze the interaction of an electromagnetic spike (one cycle) with 
a thin layer of ferroelectric medium with two equilibrium states. The 
model is the set of Maxwell equations coupled to the undamped Landau- 
Khalatnikov equation, where we do not assume slowly varying envelopes. 
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From linear scattering theory, we show that low amplitude pulses can be 
completely reflected by the medium. Large amplitude pulses can switch 
the ferroelectric. Using numerical simulations and analysis, we study this 
switching for long and short pulses, estimate the switching times and 
provide useful information for experiments. 

1 Introduction 

Recently there has been a revival of interest in the study of phenomena occurring 
during the propagation of electromagnetic waves in different media having a 
long-range order. The propagation of electromagnetic pulses (in particular, of 
optical range) as stable solitary waves, sometimes called solitons, occupies a 
special place in such investigations. As an example, we can point out to studies 
of electromagnetic solitons propagating in ferro(antiferro)magnetic media 

Based on the Landau-Lifshitz equations for magnetization and the Maxwell 
equations for the electromagnetic field, solitary waves (supersonic, as distinct 
from magnetic solitons) were found. These can be approximately described by 
the modified Korteweg-de Vries (mKdV) equation. It was shown that they are 
stable against structural perturbations of the mKdV equation caused by energy 
dissipation. Using the multiscale perturbation theory, electromagnetic solitons 
were considered in antiferromagnetic [2] and anisotropic ferromagnetic media 
[SJ. It was shown in |S] that the modulation of an electromagnetic wave in a 
fcrromagnet in an external field can be described by the nonlinear Schrodinger 
equation. This opens the way to study the modulational instability and the 
formation of electromagnetic solitons in magneto-ordered media. 

Dielectrics having a permanent polarization in the absence of an external 
electric field, called pyroelectrics, are, in many respects, similar to magneto- 
ordered media. If the phase transition between the pyroelectric and the nonpy- 
roelectric states is of second-order, such pyroelectrics are called ferroelectrics. 
The Landau phenomenological theory of phase transitions is based on the as- 
sumption of an order parameter that is zero in one phase and is nonzero in the 
other phase. A uniaxial ferroelectric provides a simple example of this theory of 
phase transitions. Due to the nonlinearity of the free energy of the ferroelectric 
relative to the order parameter, there can exist nonlinear waves of spontaneous 
polarization both in the form of solitons [7] -|10j and in the form of domain 
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walls Thin films of ferroelectrics |121 ll.ij and liquid crystals having ferro- 
electric properties ^] 1151 are of interest for practical use. This is because 
such media, having a permanent dipole moment, can be used to generate effi- 
ciently optical harmonics . Ferroelectrics can also be used to create memory 
devices and optically controlled switches El El El El ■ Many recent results 
are collected in the review on two dimensional ferroelectrics and thin polymer 
ferroelectric films represents in |19| . 

Recently [23] we investigated the traveling wave solutions of the Maxwell- 
Duffing homogeneous system describing the interaction of short electromagnetic 
pulses with a bulk ferroelectric. The description of the fast switching and ultra- 
short pulse propagation requires that the duration of these pulses be shorter than 
the relaxation time of the nonequilibrium polarization, which, for ferroelectrics, 
is equal to several nanoseconds El- In this study, we consider the interaction 
of an extremely short light pulse with a thin film of dielectric medium having a 
spontaneous polarization (a ferroelectric). Using the phenomenological Landau- 
Khalatnikov model [201 EH HU EH describing a uniaxial ferroelectric and the 
Maxwell equation for an electromagnetic wave, we consider the reflections and 
refractions of a short electromagnetic spike (i.e., a pulse without a carrier wave) 
through a ferroelectric thin film. We assumed that the film width is less than the 
spatial size of the spike, but is greater then the critical length L c below which 
no ferroelectricity exists El GUI El- According to theoretical estimations and 
E31 El L c = 0.5nm for BaTi0 3 , L c = 20nm for PbTi0 3 . For experimental 
films of two or five monolayers of polyvinylidene fluoride (PVDF) and copolymer 
P(VDF-TrFE) L c « 7nm El 

The article is organized in the following way. Section 2 presents the model 
together with a solution of the Maxwell equations for a localized ferroelectric 
medium. We apply these results to a thin film in section 3. Section 4 is devoted 
to the switching caused by large amplitude electromagnetic pulses. We conclude 
in section 5. 

2 Phenomenology of ferroelectricity 

A phenomenological description of ferroelectricity due to Landau and Khalat- 
nikov |2U1 1211 E] gives the following Lagrangian density for the interaction of 
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an electromagnetic field and a dielectric medium 

1 1 f)\~\ 

(1) 



f(f )2 -i,(P)-ipf 

2g at g c at 



where A is the vector potential, P is the polarization of the medium, c the 
speed of light, g is coupling constant, $(P) is the thermodynamic potential, the 
last term describes the coupling between A and P. The ferroelectric medium 
can exist only in some region and therefore we have introduced in /(x), 
the characteristic function of the medium, i.e., /(x) = 1 inside the medium 
and /(x) = outside. The Lagrangian approach guarantees that we take into 
account the correct couplings. 

2.1 Homogeneous case 

We first consider the homogeneous situation i.e. we assume /(x) = 1. The 
variation of the action functional yields the equations for A and P 



rotrotA 



1 d 2 A 4irdP 



c 2 dt 2 c 2 dt ' 
<9 2 P <J$(P) g OA 
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which when written in terms of the electric field E = — (l/c)dA/dt are 



1 d 2 E _ 4tt d 2 P 
d 2 P <5$(P) 



-trotE +-— = -- — , (2) 



dt 2 ' ^P ^ 
We can introduce a phcnomcnological relaxation of the medium polarization 
by adding to the right hand side of a linear damping term and obtain the 
well-known Landau-Khalatnikov equation. 

To simplify the problem we will consider an electric field polarized along the 
propagation variable x so that E = Ee x and a polarization along x P = Pe x . 
Following Landau the potential can be chosen as 

1 1 1 BP 

$(P) = $ + -aP 2 + i /3P 4 + ^D(-) 2 . (4) 

In the theory of Landau a — a^{T — T c ) depends on the temperature and 
not j3. If a > the potential is minimum for P — and this corresponds 
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to a (disordered) paraelectric phase. On the contrary if a < and > 
there are two minima located at P = ±a/ '—a/ ' (3 corresponding to two opposite 
orientations of the polarization, this is the ferroelectric phase. Then the Euler- 
Lagrange equations for E, P are 

1 d 2 E d 2 E 4tt d 2 P 



c 2 dt 2 dx 2 c 2 dt 2 

d 2 P „d 2 P . „ . „„ 3 



(5) 



D—+aP + pP J = gE 



dt 2 dx 2 

Traveling wave solutions for the system were found in |24|. they include 
a one-parameter family of solitons, a one-parameter family of kinks (domain 
wall) solutions and a one-parameter family of periodic (cnoidal) waves. In |24| 
wc found numerically that the bright solitons on a zero and nonzero polariza- 
tion background are stable while the dark solitons are stable only on a zero 
background. 

2.2 Inhomogeneous case 

We will now normalize the fields and variables as 



i' = i/VR. z' = aVR/c, A = 2 c y^a, P = J^q, (6) 

where the polarization P is normalized by the saturation value. The Lagrangian 
density 0) becomes 



l,0q, 2 q 2 1 4 9a 



2 y dt' 2 4^ 

where the primes have been dropped. The parameters are 



(7) 



7 = 2 /^ , m = a/\a\, (8) 
V l a l 

where in is the " mass" of the excitations which can be negative for ferroelectrics 
and positive for paraelectrics and 7 is the polarisability of the medium. In 
the theory of Landau-Khalatnikov, the parameter g is the susceptibility of the 
material so g = l/(4|a|) in the ferroelectric phase. From © we have 

7T 1 pi: 1 .J 

2R = V2R |Tc_T| • 
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Consider for example a crystal BaTiOz 1 then olq = 6 x 10 e K 1 , j3 = 2 x 
10~ 15 to 3 J" 1 so that 7 = 2 x 10 5 |T C - T\-\ 

The simplified Lagrangian density |2(:i| corresponding to a one-dimensional 
plane wave: 

a? ,, , ,qf q 2 q A 

£=^--f+I(x)(^-+m'L- + ^- lq a t ), (9) 

where a is the analog of vector potential. The Eulcr-Lagrange equations are 

<Ht - a xx = -jl(x)q t (10) 
qtt + rnq + q 3 = 7a* (11) 

The equations for the electric field e = —at and medium variable can then be 
obtained 

ett - e xx = -jl(x)q t t, (12) 
q tt + inq + q 3 = je, (13) 

where the coupling between the fields e and q only occurs in the medium i.e. on 
the support of I(x). 

We now consider the general scattering formalism assuming a localized elec- 
tromagnetic wave impinging on the medium from the left like is shown in Fig. 
n The general problem can only be treated numerically so we simplify it and 
consider two limiting cases, an array of thin films and a single thin film. 

2.3 Scattering of an electromagnetic wave by a ferroelec- 
tric slab 

We assume that the electromagnetic wave is incident from the left x < on a 
medium whose position is given by the indicator function I(x). Then the wave 
equation reads 

ett ~e X x = g{x,t) (14) 
g(x,t) = -jq tt I(x), (15) 

where the boundary conditions are 

e(t,x — ±00) = 0, et(t,x = ±00) = 0, (16) 
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Figure 1: Schematic of an incident electromagnetic pulse on a thin ferroelectric 
film panel a (left). The b and c panels show respectively the case of switching 
and non switching of the polarization. 

and the initial conditions 

d 

e(t = 0,x) = e (x), e t (t = 0,x) = e x (x) = -— e {x). (17) 

We suppose that the initial pulse eo(x) is located far to the left of the medium. 
The equation for the field is linear so to solve we introduce the Fourier transform 
of e and similarly for g 

+00 +00 

/dk f 
— — e(k, t) exp{ikx), e(k,t) — / e(fc, t) exp(— ikx) dx, (18) 

— 00 —00 
to obtain the initial value problem 

e tt + k 2 e = g(k,t), (19) 

with e(k,t = 0) = eo(fc), et(k,t = 0) = ei(fc). 

Following the general approach 2Ej we write e = e° + e 1 where e° solves the 
homogeneous equation 

e° + k 2 e° = 0, (20) 

with the initial conditions e°(k,t = 0) = eo(fc), e^(k,t = 0) = ei(fc). and e 1 
solves the inhomogeneous equation 

si , 1.2-1 



k z e l =g(k,t), (21) 
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with e x (k,t = 0) = 0, e\{k, t = 0) = 0. 
The general solution of 1|2U|) is 



1 



= 2 ^o(fc) + ^i( fc )J exp(-ikt) + - ^e (fc) - -j^x{k) J cxp(+ifci), 

Using for example a Green's function approach one easily sees that the solution 
of the problem l|21[l is 

t 



sin k(t — t) 



g(k,T)dr, 



The solution of the homogeneous problem can be transformed into the fol- 
lowing expression 

x-t 

e °0M) = ^ e (x - t) + e (x + t) + J ei(y)dy = e (x - 1) 

x+t 

the first part of which is D'Alembert's formula. The inhomogeneous solution 
can be rewritten as 

+oo t 

e 1 (x,t) = ^ J J g(y,r) [6(x - y - 1 + r) - 6{x - y + 1 - r)] drdy, 

-oo 

z 

where we use the step-function 6{z) = J S(x)dx. So we have the general 

— oo 

solution of the scattering problem under consideration. 



-t-oc t 



e(x, t)=e {x-t) + - I I g{y, r) [6[x - y - t + r) - 6(x - y + t - r)] drdy 



-oo 



(22) 

As an example we give the result for an array of thin films, useful for appli- 
cations and theory. The indicator function is 



JY 



I(x) = y^ j 5(x- Xi), 



(23) 



and the final result is 



e(x, t) = e (x — t) + ~^2 I ' 1t(xi,r) [5(t - t + x - x t ) + 5(t - t - x + x t )] dr, 
i=i 

(24) 
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where the evolution of the state of every film at point x = x% is defined by the 
equation 

q u {xi, t) + mq(xi, t) + q(xi, t) 3 = ~/e(x = x iy t). 

This problem is still complicated to analyze because of the delays and advances 
in equation l|24l) . We therefore simplify drastically the situation by considering 
a single thin film only. 

3 A single thin film 

There is the simplest not trivial case when the nonlinear medium is represented 
by a single thin film of anharmonical oscillators. Here we can calculate the 
integral in expression (|24[1 exactly. The electric field is defined by the following 
expression 

* 

e(x,t) = eo(x-t) + lfq t (0,T)[S(T-t + x) + S(T-t-x)]dT. (25) 

o 

At the point where the film is placed [x = 0) we have 

e(0, t) = e (-t) + | j qt(0, r) [S(t - t) + S(t - t)] dr 
o 

= eo(-t)-|gt(0,t). 

Substituting this expression into the equation of the oscillator results in the 
following equation 

q tt + mq + q 3 = je (-t) - —q t . (26) 

This model represents the evolution of the nonlinear Duffing oscillator with 
both damping and forcing. 

3.1 Linear considerations 

In the homogeneous case g = 1 we can compute the dispersion relation for the 
system (|1U|I by assuming a = aoe l ( kx ~ ut \ q — q e llykx ~ wt ^ and obtain 

k = 2 +uj 

The dispersion relation u>(k) presents the well known polaritonic gap. 
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When the medium is localized like for the thin film case, the picture changes 
and we need to make a scattering experiment to understand the linear behavior 
of the device. The equations (|12|) become for the paraelectric case 

e t t~e xx = -jq tt S(x), (27) 
qtt+q = je(x = 0,t). (28) 

We compute the reflection and transmission coefficients of an incident linear 
wave on such a medium. For that we introduce the harmonic dependence 
e(x,i) = e luJt E(x), q(t) — Qe %U3t . Plugging these into the previous system 
of equations yields the Schroedinger equation with delta function potential 

,2 

E = 0. (29) 



E 



XX 



" 2 + T^ 



For the scattering we assume a wave incident from the left x < E — e^ lkx + 
Re lkx and a transmitted wave x > E = Te~ lkx . The continuity of E at 
x = and the jump of the derivative [JSJL = — 7^ 2 i?(0)/(l — w 2 ) give the 
transmission and reflection coefficients 

2i(l 

2i(l-w 2 )-7Cj' 



T '—, (30) 



R = T-l = — 1UJ , (31) 

2^(l — u z ) — juj 

where we used the dispersion relation k = u>. 

Several remarks can be made. First transparency R = is obtained only for 

uj = and total reflection occurs for cu 2 — 1 as expected [35] ■ We also have two 

bound states corresponding to the poles of R and T 



which are located in the lower half complex plane. 

In the ferroelectric case the medium will oscillate around one of the equilibria 
q = ±1 so that the linearized equations are 

e tt -e xx = -jq tt 5(x), (33) 
q tt + 2q = je(x = 0,t). (34) 

The reflection and transmission coefficients become 

R = T — 1 = 7L 1 . (35) 

2i(2-w 2 ) -juj K ' 
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Figure 2: Modulus of the reflection coefficient \R(uj,j 2 )\ 2 . 
In this case the poles are given by 

^ J -T ± f I f,- < 36 > 

In Fig. |21 we show the reflection coefficient \R(u>, 7 2 )| 2 - It is close to 1 for 
the resonant frequency uj — 1 (resp. \/2) in the paraelectric (resp. ferroelec- 
tric) case. The width of the resonance is the same for both cases and increases 
when the coupling 7 grows. The coefficient 7 is the coupling between the elec- 
tromagnetic field and the medium. It is appears as a damping of the medium 
polarization in l|26|) . The mechanism of relaxation is radiative. 



4 Effect of anharmonicity: switching 

When the amplitude of the incident pulse is large enough the anharmonic term 
in H26|) needs to be taken into account. If the system is paraelectric (to > 0) then 
it gets kicked out of the equilibrium state q = and relaxes back to it. More 
interesting is the ferroelectric system (to < 0) which has two stable equilibria 
q = ±1 and one unstable q = so that switching between them is possible. We 
will solve 

qtt ~ q + q 3 = je {-t) - — g t , (37) 

for an incoming electromagnetic pulse eo(— t) assuming the ferroelectric medium 
is in its equilibrium position i.e. with the initial condition qt(t — —00) = 
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0, q(t = — oo) = q = — 1. For a harmonic perturbation, chaotic behavior can 
occur, see for example |27) . Here we have a force of finite duration so we expect 
transient chaos which will cause irregular switching. 

From the experimental point of view it easier to work with pulses whose 
profiles are gaussian or plateau-like form. We chose the plateau-like form 

e (-t) = A Q ( tanh ~* - ^ - tanh — - - ) . (38) 

V */ */ / 

The initial polarization of the ferroelectric medium is defined by the parameter 
q(t = -oo) = q. 

4.1 Short electromagnetic pulse 

For a short electromagnetic pulse the equilibrium positions can be considered 
as fixed. After the interaction, the system evolves as if free. Figure shows 
the evolution of the medium polarization and the corresponding phase plane 
under the action of an ultrashort electric field. The amplitude is not enough for 
switching and the system relaxes to q = — 1 following the linearized behavior 
exp(iu>t) where lu is given by (|36f) . If the amplitude is increased as in Fig. [3] the 
system switches to the fixed point 5=1 following again the linearized behavior 

4.2 Slowdown of switching 

Now we consider the special cases where the incoming pulse brings the system 
on a trajectory that goes near the unstable fixed point q — 0. Fig. shows 
such a case below threshold. The left panel of Fig. shows the slowing down of 
q(t) around q — because there dq/dt sa 0. As seen from the plots, the system 
remains "frozen" near the unstable equilibrium state q = for a certain time 
Tdei- In Fig. |B|we present the results for a slightly duration where the system 
has switched. 

To investigate this slowing down, we introduced a gaussian initial electric 
field 

e (-i) = A Q exp[-t 2 /a 2 }. 

In Fig. [7| we plot a typical evolution q(t) indicating Tdei on the left panel. On 
the right panel we give log(Td e i vs log(P) where P is the total power in the 
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Figure 3: Plot of q(t) (left panel) and phase portrait (q,dq/dt) (right panel) 
for a incident pulse below threshold so that the medium does not switch. The 
parameters are A = 0.2, t\ — —20, t 2 — —22, tj = 0.5, 7=1. 




Figure 4: Same as fig. |3]but with a larger amplitude Ao = 0.4 so that switching 
occurs. 
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Figure 7: Plot of q(t) (left panel) showing the slowing down as the system 
approaches the unstable fixed point q = 0. The right panel shows the duration 
of the plateau Tdei vs the power P of the incident pulse in log-log scale 



electromagnetic pulse 

P = J +0 ° dte 2 {t) = Ala^. 

We present three different amplitudes Aq = 1, 0.75 and 0.5 and vary P by 
varying the width a. Fig. shows that there are critical values of power where 
the time delay becomes very large. There does not seem to be a single critical 
index for the description of the singular behavior of the time delay. 

4.3 Free evolution of the oscillator 

We consider the limit when the external electromagnetic field is so short that 
it just gives an impulse to the oscillator which then evolves as free. Then 
eo(t) — AoS(t). Plugging this into (|37|l . integrating in a small interval around 
t — and assuming continuity of q we obtain qt\t=o+ = 7^0- The free evolution 
of the oscillator is then governed by the equations 

dq/dt = p, 

dp/dt = q- q 3 - 0.5j 2 p, 
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Figure 8: Phase portrait (q, q) showing free evolution of the system for a delta 
function like electromagnetic pulse so that q(t — 0) = —1, dq/dt(t = 0) = 
—0.75. The long dash corresponds to 7 = 0.2 and the short dash to 7 = 0.5. 
The continuous curve corresponds to 7 = 0. 

with the initial conditions q(t = 0) = — 1, q t (t = 0) = 7A0. The phase 
trajectories are solution of the equation 

p (dp/dq) = q — q — 0.5j 2 p. 

The phase planes corresponding to this free evolution of the Duffing oscillator for 
7 = 0.5 and 0.2 are presented in Fig. [^together with the nondamped situation. 
For 7 = 0.5 the damping is strong and prevents switching while for 7 = 0.2 
the system can escape to the other equilibrium and slowly converge to it. The 
pictures on Fig. [SJshow that dissipation leads to damping of oscillation around 
equilibrium position and supports the switching from one to another equilibrium 
state. This dissipation results from the radiation of the electromagnetic waves 
out of the thin film. 

4.4 Long electromagnetic pulse 

When the electromagnetic pulse has the form of a plateau with a sharp front 
and a sharp tail, transient steady states of the ferroelectric are created. We can 
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find them by considering the static solutions of equation (|37|) : 

q-q 3 + 7 e = 0. (39) 

Let eo be positive. If < 760 < 2/(3\/3) there are three roots, corresponding 
each to a fixed point, two stable and one unstable. These can be calculated by 
perturbation when jeo << 2/(3v3), we have 

<Zo(eo) ~ -7^0, 9i,2(eo) « ±1 + 7 e o/2. (40) 

When jeo = 2/ (3v3) the unstable and left stable equilibrium points merge 
together and we have only one stable point 
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(e ) = 2/V3. (41) 



In the limit of very big amplitude of electromagnetic pulse, when ^cq >> 
2/(3-^/3), there is only one fixed point which is stable and is defined by the 
approximate formula 



<72(e ) ~ (7 e o) 1/3 



1 I \-2/3 

2 (7eo) 



(42) 



The numerical simulation of switching under the influence of long plateau- 
like pulse demonstrates the damping nutations near the stable points (71,2(60). 
The kinetics of switching can be considered by using the linearized equation 
IpTfjl near the stable fixed points: 

7 2 

Sqtt + Y 5 1t + [ 3 9i,2(eo) - 1]Sq = 0, 

where 5q = q — 91,2(60)- The associated characteristic equation give the decre- 
ment and frequency of nutations 



r = 7 2 /4, fio = ^/3 g t 2 (eo)-l-7 4 /16. (43) 

The decrement is independent of the initial pulse amplitude while the frequency 
of nutation depends on it. 

The following evolution of the ferroelectric polarization depends on the area 
of the external field. In this section we consider long pulses that create transient 
equilibria in the ferroelectric. Fig. 03 shows in the left panel such a long pulse, 
in the middle panel the response q(t) of the medium and in the right panel the 
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Figure 9: Plot of a long incident pulse eo(x) (left panel), q(t) (middle panel) 
and phase portrait (g,dq/dt) (right panel). The parameters are the same as in 
Fig. EJexcept A Q = 0.1, t x = -70, t 2 = -20. 




20 40 60 80 100 120 -1 1 2 

t q 

Figure 10: Plot of q(t) (left panel) and phase portrait (q, dq/dt) for an incoming 
long pulse similar to the one in Fig. except that A = 0.5 
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Figure 11: Parameter plane (amplitude, duration) of incident pulse detailing the 
final state of the film switched or non switched. The + symbols correspond to 
no switching and the x to switching. 



associated phase plane. The ferroelectric is moved into a transient equilibrium 
and returns to its previous stable polarization state q = — 1. Here we chose 
Aq = 0.1 so that 7eo = 0.2 in the plateau region. The system goes to the 
transient steady state q~\ w —0.88 value that is in good agreement with (|40() . 
which gives —0.9. The decrement Tq ~ 0.25 and nutation frequency f2o ~ 2 is of 
the motion to this transient fixed point is predicted correctly by the estimates 
Ij43(l . When the system returns to its natural fixed point, the estimates are 
again correct and given by l|36[) . Notice that the nutation frequency around the 
transient fixed point is twice as big as the one around the natural fixed point. 

In Fig. HOI we observe the same phenomenon except that the transient equi- 
librium is close to the other stable polarization state q — +1 so that when the 
field returns to zero, the system relaxes to that state. Here Aq — 0.5,7 = 1 so 
that in the plateau region eo = 1. We are in the region above the critical 6q so 
that there is only one fixed point. The estimate 142|) gives q « 1.5 which is in 
excellent agreement with the numerical value 1.4. 
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5 Concluding remarks 



We considered the simplest model of interaction of a short electromagnetic pulse 
with a thin ferroelectric medium where the polarisation can be considered as 
uniform. The duration of the electromagnetic pulse is much shorter than the 
relaxation time of the medium so that only radiative decay occurs. 

The linear scattering formalism predicts that low amplitude pulses can be 
completely reflected by the medium, in a frequency band that grows with the 
coupling 7. On the contrary strong electromagnetic fields can switch the medium 
from one state to another. We have studied such switching phenomena for both 
short and long pulses. For the latter we characterized the transient states they 
create. 

We define the switching time t s as the time interval for the system to go from 
one fixed point to the other. From the Duffing system the typical (normalized) 
damping time is 7 2 /4. If j 2 > 4 switching occurs during the front of the pulse, 
otherwise the switching time is of the order of 7 2 /4 because the system circles 
around the fixed point before reaching it. In any case the system will always 
switch in a time smaller or of the order of 7 2 /4. In physical units this is about 

t s =A\a Q \~^\T c -T\~V\ 

In some cases, the switching time can be longer, in particular if the field drives 
the ferroelectric near the unstable fixed point causing considerable slowing down. 
Switching is also irregular as shown by Fig. 1111 which gives the events in the 
plane (duration, amplitude) of the incoming pulse. There one sees that a thresh- 
old amplitude is needed for switching. Above that the system switches or not de- 
pending on the pulse duration. For a given duration there seems to be windows 
where the system switches. All this information can be used by experimentalists 
to estimate parameters. Finally we believe the model due to its simplicity and 
generality can be transposed to other electromagnetic systems. 
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